Berry phase in open quantum systems: a quantum Langevin equation approach 
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The evolution of a two level system with a slowly varying Hamiltonian, modeled as a spin 1/2 
r**"* ' in a slowly varying magnetic field, and interacting with a quantum environment, modeled as a 

bath of harmonic oscillators is analyzed using a quantum Langevin approach. This allows to easily 
obtain the dissipation time and the correction to the Berry phase in the case of an adiabatic cyclic 
evolution. 
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PACS numbers: 

03.65.Vf Phases: geometric; dynamic or topological 
03.65.Yz Decoherence; open systems; quantum statistical methods 

03.67.Pp Quantum error correction and other methods for protection against decoherence 

I. INTRODUCTION 

In his seminal work [l[ Berry showed the appearance of a purely geometrical phase factor associated to the non 
degenerate eigenstates of a Hamiltonian undergoing a cyclic adiabatic evolution. Since then much work has been 
done to generalize such concept to non cyclic, non degenerate or non adiabatic evolution [2J • The renewed interest for 
geometric phases in a quantum computation scenario @, 0, H, ||| 0] is due to their supposed intrinsic fault tolerance. 
Such hypothesis has been analyzed in for the case of a spin in a slowly varying magnetic field with small classical 
random fluctuations. There it was shown that for small fluctuations, i.e. to first order in the perturbation, and in the 
adiabatic limit the main source of decoherence are dynamical fluctuations. Similar conclusions have been reached for 
quantum noise in [9( by means of a quantum trajectories approach. In the present paper we will analyze the problem 
using a quantum Langevin equation approach. Our system consists of a pseudospin interacting with a quantized 
bosonic field. The spin free Hamiltonian is assumed to undergo a slow cyclic evolution. The geometric phase appears 
^J^" in a natural way in terms of the so called adiabatic Hamiltonian [lOj . Once such Hamiltonian is introduced the 
Heisenberg equations of motion are derived and from them the quantum Langevin equations. This approach allowed 
• i-h , us to analyze the effects of the quantum fluctuations on both the decay constants and on the overall phase acquired by 
the spin energy eigenstates in their cyclic evolution. The same problem has been addressed by some recent papers [l2| 
with the use of an elaborated perturbative technique. We will show how our approach allows to obtain in a simpler 
and more straightforward way the corrections to the the Berry phase found in [121 ] ■ Furthermore we will provide a 
transparent physical picture of the results obtained. The effect of noise on geometric phases in different scenarios 
from the one described above has been studied in 
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II. THE ADIABATIC HAMILTONIAN 



The system we consider consists of a pseudospin in a slowly varying static magnetic field interacting with an 
environment modeled as a bath of harmonic oscillators. The overall system Hamiltonian is assumed to be of the 
standard form 

H = ■ a + ^tu k ala k +^g k a z (a k + a{^j (1) 
k k 

where 3 = (a x , o~ y , a z ) are Pauli operators and B(i) = £?o(sin $ cos tp, sin$ sintp, cos $) is a three dimensional vector, 
which we assume to be time dependent and a k (a k ) are bosonic creation (annihilation) operators for mode k (we have 
set h = 1). 
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The first step to obtain the Heisenberg equation of motion for the spin and bath operators is the introduction of 
the so called adiabatic Hamiltonian [IcJ i.e. of the Hamiltonian whose eigenstates, in the absence of interaction with 
the environment, after a cyclic evolution acquire the dynamical and geometrical phase predicted by Berry. Let us 
first rewrite the free spin Hamiltonian in the form 



H S ee -B • a = -f (| U (t))(U (t)\ - | U (t))(i n (t)\) (2) 

where | t„ (t)) and | J.„ (t)) are the eigenstates of Hs at time t i.e. the eigenstates of the operator a ■ n, where 
n ee (sin d cosy, sin ■& sin ip, cos I?) is a unit vector pointing in the instantaneous B direction. Let us then define the 
following time dependent unitary operator: 



U(t) = \U(0))(U(t)\ + \i n (0))(i n (t)\ (3) 

In the absence of any coupling with the environment the time evolution of the state vector \ip(t)) = U(t)\ip(t)) is 
generated by the Hamiltonian 



H S = U(t)H s U\t) - iU(t)j t U\t) (4) 

A time dependent B will in general induce transitions between the instantaneous energy eigenstates | ] n (t)) and 
i n (t)). However if the direction of B changes slowly enough in time we can neglect such transitions, which in (|3|) 
are described by the term j(| i (t)). This is nothing but the standard adiabatic approximation [llj], valid as 

long as i(t {t)\f t \ I (*)) < B. On the other hand i{] (t)\j- t \ ] {£)) and i{[ (t)\-^\ | (t)), which are nothing but the so 
called Berry connection, are responsible of the appearance of a non vanishing geometric phase and therefore must not 
be neglected. The adiabatic Hamiltonian is therefore 

Hf= (^-*'<Tn(*)|^| T»(t)>) I Tn(0)>(t n (0)|- (j-+i(Ut)\j t \ U(*)>) I U(0))(1„(0)| (5) 

Note that when B undergoes a cyclic evolution the eigenstates of ([5]) correctly acquire the dynamical plus the geo- 
metrical phase predicted by Berry. 

In the basis of the eigenstates of the a z operator we have 

|Tn) = e-^ 2 cos^|t z }+e^/ 2 sin^|I z } 

| U = e-^/ 2 sin|| U -e^ 2 cos 1 1 | z ) (6) 

from which it follows that 



i(U{t)\^\ T„(i)) = -*Un(*)l^l InW) = 4 C ° S,? ^ 



The adiabatic Hamiltonian takes therefore the form 



^= *°~7 S ^W * (8) 

For the sake of simplicity, and with no loss of generality, we will consider the case most discussed in literature in 
which B(i) precesses slowly with angular velocity = 2tt/T i.e. (p(t) = fit, (p(0) = 0,^(i) = $(0). Furthermore we 
will rotate our axis so that [cosi?(0)cr z + sini?(0)CTj;] — > cr z In this case the adiabatic Hamiltonian takes the following 
simple form 



frad _ u o 

H s - Y z 



(9) 
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where uiq = Bq — ficos$(0). In order to simplify the comparison between our results and the existing literature we 
would like to point out that our approach is different from, although at the end equivalent to, the one in which a 
rotating frame is introduced (see. e.g. [II]). In geometric terms this implies that a field of amplitude along the z 
direction is added to B, i.e. in the rotating frame the effective magnetic field changes both in direction and length. 
This is not the case in our adiabatic Hamiltonian, in which the magnetic field changes only in length. However when 
the overall accumulated phase is calculated to lowest order in Q their model coincides with ours. In other words 
while we assume from the very beginnings the adiabatic limit in [l2| the non adiabatic contributions are discarded a 
posteriori. 

To conclude this section we write our full, time independent, adiabatic spin - boson Hamiltonian 



H ad = ^-<7 Z +'^2ui(k)ala k + ~^9k (&k + a\) {cr z cost) - o x sin??) 



(10) 



Note that in (|10[) and from now on the spin operators are in the adiabatic frame. 



III. THE LANGEVIN EQUATIONS 



The quantum Langevin equation of motion for a two level systems coupled with a harmonic bath is well known in 
literature Here, for the sake of clarity, we will sketch the main steps of its derivation. The starting point are 

the Heisenberg equations of motion (O = i[H, O]) for the spin Pauli operators o>, er+ = (a x + ia y )/2 and the bath 
operators a k , a k which are 



& z = 2i sin ■& g k (cr + - <r_ ) (a k + aj,,) 

k 

<r+ = tcjQ(T + + g k (a k + a];) (sin'i9(T z + 2 cost?q+) 

k 

ak = -Wk&k — igk (cosi?cr z - sin#(cr+ + cr_)) 



(11) 
(12) 
(13) 



which can be cast in the following integral form 



<J z (t) = ( 7 z (0) + 2^sintf^ 5fc / dt' [a+(t') - a_(t')] fo*(f) + 4(*') 



o+(t) = a + (0)e^ at +iY,9k / dt'e 



O iuo(t-t ) 



a k (t') + at(i') [sintfcr z (i') + 2costfcr + (t')] 



a k (t) = a k (0)e 



—iuJkt 



ig k I dtfe-*"^*-* } {cosi?CT i (i')sinT?[<T+(i / ) +<7_(t')]} 
'o 



(14) 
(15) 
(16) 



A standard assumption in the derivation of a quantum Langevin equation is that the timescalc of the decay processes 
is much slower than the free evolution. In other words the weak coupling with the bath degrees of freedom modifies 
the spin dynamics on timescales which are much longer that uJq 1 . In the above integrals we can therefore put 



Cz(t ) = cr z (t) 

o+(t) = e- ia ">(*-*'V + (t) 



Ofc(t 



iu) k (t-t )- 



a k (t) 



Eqs. dligilBJ then become 

cr z = <J Z (0) + sini?^g fe |a+a fc C*(wo - Uk) + v+a{Ci UJ o +^fe) - <7-a fc C(^o +^fe) - (T_a[C(^o -Uk) 



(17) 
(18) 
(19) 

(20) 



a, = e tuJot ' 



cr+(0) + i sin d'^2g k a z a k ((u>o + u) k ) + %C(^o ~ Wk) + 2icosi?^g fc CT+ a k C,{uj k ) + a[(*(wfc) (21) 



ak 



-iuj k t 



a* (0) -ig k [cos $o z C (wfc) - sin ■& (a + (* (w + ) + <*- C (wo - u k ) ) ] 



(22) 
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where 

/•*.'< i 

C{x) = iim / e rxt dt =P- + irS(x) (23) 

P denotes principal part and S(x) is the Dirac delta function. The integration limit t — > oo is justified on the ground 
that we are interested at times t 3> 0Jq . 

Inserting the (f2"0"l - 12"2"| into the (JTT]- [T3J) and taking care of a consistent choice of all operator products [H|,[l6[, one 
obtains the desired equations of motion for the spin operator in the rotating frame, averaged over the environment 
degrees of freedom 



= -2sin 2 ^(7_Lcr z +7^„ ac l ) - 2sin2?97|| + cr_) 



(24) 



dcr+ 



[2(i£-7_L„ac)l - («A + 7i)fT z ] +sin 2 i9{(iA-7±)cr + - + 7_l)cx-} - 4cos 2 •d^ap.h) 



where, transforming the sums into integrals 

7_L = 7T / duj k p(ui k )gl (2n k + 1) S(uj - ui k ) 
Jo 

/>LU C 

7|| = 7r / duj k p(u k )gl(2n k + l)6(uj k ) 



A 



dui k p{ui k )gl{2n k + 1) 



P 



P 



dui k p(ui k )gl 



2P 



LU —U)k ^0 + u k 

p p 



LO k LU - LU k LU + LO k 



(26) 
(27) 
(28) 
(29) 



In the above equations p{uj k ) is the density of modes at frequency Lu k , n k is the mean number of photon in field mode 
k, 1 is the identity in C 2 and jj_ vac is j± for n k = . 



IV. DISSIPATION AND ENERGY SHIFTS 



The above equations allow us to clearly identify the effects of the adiabatic evolution on the physical quantities 
which characterize the spin dynamics, namely the decay constants and the energy shift. First of all let us consider 
the decay constant 711 which describes the decoherence mechanism due to fluctuations "parallel" to the instantaneous 
direction of B. As expected it is not modified by the adiabatic change of such direction. Furthermore the value of 711 
depends on the density of field modes at zero frequency which, in most situations of physical interest is equal to zero. 

The decay constant 7^ describes the dissipation mechanism due to the exchange of energy between system and 
bath and depends on the density of modes at the resonance frequency uo = ujo ■ If we assume that the density of modes 
is a slowly varying function of to near resonance, we can safely assume for very small f2, i.e. in the adiabatic limit, 
J " du} k p(u! k )g1 (2n k + 1) 8(u> — u> ) ~ J Q C duj k p(uj k )g1 (2n k + 1) 5(u> — B ). This confirms that the timescale of dipole 
decay is not modified by the adiabatic evolution, a result which has been obtained with different techniques, from 
classical stochastic noise Q, to quantum jump Q. We should point out that in order to observe the geometric phase 
we must have 



7^ <C ft <C ujq 



(30) 



Let us consider now the change in the energy shift A. In the adiabatic limit we must consider terms up to order 
0(tt) and therefore 



A = / dio k p(ijj k )g k (2n k + 1) 
'0 



P P \ d 

Bq — oj Bo + Lu I oojq 



p 



p 



w A + 6X 

The quantity sin 2 $Aq is nothing but the Lamb Shift [l5|, [lj| , while 



(31) 
(32) 
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SX = ftcostf / du> k p(ui k )gl(2n k + 1) 

Jo 

gives information on the the effect of the quantum fluctuations on the geometric phase. This correction coincides 
with the results obtained by [r| with an elaborated perturbation technique. The observable overall phase difference 
between the two energy eigenstates at the end of their cyclic evolution, i.e. at time T = 2tt£1~ 1 , will be 



1 



1 



{B -uY T (B + uf 



(33) 



$(T) = <$> D + $ G 



(34) 



where the dynamical phase $£> 



$ D = 



B a + sin til du; k p(u; k )g k (2n k + 1) 



P 



P 



Bu 



Bq+oj 



T 



is simply due to the renormalized energy splitting, while the geometric phase $g is 



$ G = 2ir cos d i 1 



sin 



i? / du>kp(uk)gl(2n k + 1) 



1 



(S 



(S 



(35) 



(36) 



The expression (|36[) is amenable to a straightforward intuitive geometric interpretation. It is a well known fact that 
for a spin 1/2 the Berry phase is is equal to the solid angle spanned by the time varying magnetic field B on a unit 
sphere centered around degeneracy. As opposite energy eigenstates acquire opposite geometric phases the overall 
phase difference between them will be, for a slowly precessing field at an angle equal to $ Berry = 27rcosi9. In the 
presence of a weak coupling with the bosonic bath however each energy eigenstate will undergo virtual transitions, 
responsible for the Lamb Shift, with a probability 



Prob v t = sin d 



du> k p{uj k )gl(2n k + 1) 



1 



1 



(B 



(B 



(37) 



During such transition the spin state parallel (antiparallcl) to the direction of the field B 'jumps' to the antiparallcl 
(parallel) spin state, acquiring an opposite geometric phase. The overall geometric phase difference between the 
energy eigenstates will be therefore decreased by an amount proportional to Prob vtl as shown in (|36p . Notice that 
the correction to the Berry phase is of order 0(g 2 ). In where the effects of classical noise were considered, no 
analogous correction was obtained because only contributions to first order in the fluctuating field were taken into 
account. 



V. CONCLUSIONS 



In this paper we have shown how the corrections to the Berry phase and the decay constants for a spin 1/2 
undergoing an adiabatic cyclic evolution can be obtained in terms of quantum Langevin equation once the adiabatic 
Hamiltonian is introduced. We have confirmed that the main source of decoherence is due only to the dynamical 
fluctuations, a result which has been obtained with different techniques, from classical stochastic noise @, to quantum 
jump Q and which emerges in a straightforward way in our approach. The Heisenberg equations of motion give also 
the correction to the geometric part of the overall phase difference between the energy eigenstates at the end of the 
cyclic evolution due to the coupling with the bath. We have also shown how such corrections are amenable to a 
straightforward geometrical interpretation. 
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